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Tunneling electron induced rotation of a copper phthalocyanine molecule on Cu(lll) 
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The rates of a hindered molecular rotation induced by tunneling electrons are evaluated using 
scattering theory within the sudden approximation. Our approach explains the excitation of copper 
phthalocyanine molecules (CuPc) on Cu(lll) as revealed in a recent measurement of telegraph 
noise in a scanning tunneling microscopy (STM) experiment [Schaffert et at, Nat. Mat. 12, 223 
(2013)]. A complete explanation of the experimental data is performed by computing the geometry 
of the adsorbed system, its electronic structure and the energy transfer between tunneling electrons 
and the molecule's rotational degree of freedom. The results unambiguously show that tunneling 
electrons induce a frustrated rotation of the molecule. In addition, the theory determines the spatial 
distribution of the frustrated rotation excitation, confirming the striking dominance of two out of 
four molecular lobes in the observed excitation process. This lobe selectivity is attributed to the 
different hybridizations with the underlying substrate. 

PACS numbers: 68.37.Ef,72.10.-d,79.20.Rf,74.55.+v 

Keywords: inelastic tunneling, rotational excitation, molecular adsorption, telegraph noise, molecular ma- 
nipulation 
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I. INTRODUCTION 

Tunneling electrons permit to induce changes in atomic 
or molecular adsorbates at surfaces with great con- 
troPHUl. This breakthrough in fundamental science and 
in nanotechnology has stirred a lot of attention as it be- 
came possible to manipulate adsorbates and to induce re- 
actions on the atomic scale. On the one hand, these tun- 
neling experi ments provide perfect toy systems to learn 
chemical rulealSIill^ on the other hand, promising bottom- 
up techniques for creating nanometer scaled devices be- 
come feasible^. 

The first atomic manipulation experimentsHS led to 
an important theoretical effort to unravel the mecha- 
nisms behind atomic motion induced by tunneling elec- 
trons. The task was complex due to the involvement 
of several length scales and the large number of de- 
grees of freedom. Most treatments have favored master- 
equation approaches where the atomic and electronic de- 
grees of freedom are perfectly separated and electronic 
transitions are incorporated only through atomic exci- 
tation and de-excitation ratea^SlEi] These rates are ob- 
tained through full quantum mechanical calculations, us- 
ing electron-atom matrix elements fitted to or extracted 
from deformation-potential types of calculations in a 
Golden-rule treatmenpS. The atomic evolution is sup- 
posed to be that of a truncated harmonic potential, 
whe re th e truncation indicates the rupture of an atomic 
boncpESl. Only a few works have treated different atomic 
potential energy surfaces (PES) other than the truncated 
harmonic. Avouris et ar^ used atomic wave-packet 
propagations to evaluate the atomic dynamics in an ex- 
cited PES after electron tunneling. More recently, the an- 



harmonicity of the PES as revealed by density functional 
theory (DFT) was used to explain the ele ctron -induced 
motion of ammonia molecules on Cu(100pI2ll, Despite 
the simplicity of these models, much insight was gained 
for single-atom and single-molecule dynamics. 

Among the different evolutions of molecular dynamics 
on surfaces, rotations have proven to be complex'^ ^^. The 
origin of this difficulty can be traced back to the excita- 
tion mechanism itself. Instead of inducing a deformation 
of some localized bonds due to a brief charged state of 
the molecule, rotation implies some type of partial angu- 
lar momentum transfer. Indeed, rotational excitation is 

close r to mag netic excitations than to vibrational excita- 
tions.ESISiM 

By studying the telegraph noise observed in the tun- 
neling current, Schaffert and coworkerJISl have recently 
reported on the frustrated rotation of single copper ph- 
thalocyanine (CuPc) molecules on Cu(lll), excited by 
the tunneling current. Here, we present a theoretical 
account of this rotational excitation process. The adia- 
batic PES for the rotation of CuPc on Cu(lll) is com- 
puted within DFT. We characterize the molecular ad- 
sorption as well as its electronic structure and calculate 
scanning tunneling microscope (STM) constant current 
images. These predictions are compared with the experi- 
mental observations^^ and are used to explain them (Sec- 
tion In]). The electron-induced excitation is analyzed after 
the DFT structure calculations. A scattering theory ac- 
count of the rotational excitation is presented. Using the 
customary approximation of the Tcrsoff-Hamann treat- 
ment of STM images'^ we obtain a DFT-based spatially 
resolved description of the rotational excitation by tun- 
neling electrons (Section |III[ ). Finally, the insight gained 
by our theoretical treatment is discussed and the article 



ends with some concluding remarks. 



II. STUDIED SYSTEM 

Our previous experimental and theoretical worlPS 
shows that the tunneling current induces a frustrated 
rotational motion of CuPc molecules on Cu(lll). An 
in-plane molecular axis changes by '^ 7° back and forth 
around the central Cu ion, which remains located at a 
highly symmetric bridge position on the surface. The ro- 
tation corresponds to transitions between the minimum 
energy conformation aligned with the [—110] direction 
of the surface and two metastable local minima for the 
clockwise and anti-clockwise rotations. 

A brief summary on the experimental findings reads 
as follows. When CuPc (four- fold symmetric in the gas 
phase) is adsorbed on the six-fold symmetric Cu(lll) 
surface, the resulting STM images appear two-fold sym- 
metric. One pair of opposing benzopyrrole rings (lobes) 
appears pronounced as compared to the remaining two 
lobes. In the submonolayer regime, CuPc adsorbs indi- 
vidually, the formation of dimers or clusters is rather un- 
likely. In addition, the molecules appear partially fuzzy. 
The blurring of the molecular images is attributed to 
a molecular motion induced by the tunneling electrons. 
This is due to switching transitions between two discrete 
levels that correspond to two different tunneling currents 
for a fixed tip position. The switching events are ran- 
dom in time, with constant probabilities .f^^ Therefore, 
the definition of random telegraph noise (RTN) is ful- 
filled. The high-current level corresponds to nearly twice 
the low-current level. The switching frequency scales lin- 
early with the tunneling current, hence, a one-electron 
driven process is observed. 

A special electronic setup was designed to analyze the 
RTN signal from the tunneling current in real-time, dur- 
ing the ongoing STM experiments. For a full character- 
ization of a telegraph signal, the three quantities rate, 
amplitude and duty cycle have to be measured. Details 
on this new technique, scanning noise microscopy (SNM), 
are discussed in Ref. [55]. The RTN characteristics can 
be obtained in spectroscopy measurements with fixed po- 
sition of the tunneling tip, but they can also be mapped 
to create noise images exhibiting e.g. spatially resolved 
excitation rates with the same resolution as the STM 
topography. The study of CuPc on Cu(lll) yielded ex- 
citation maps clearly highlighting a lobe selectivity of 
the rotation excitation mechanism. Two out of the four 
molecular lobes appear noisy in STM. In SNM only those 
two lobes are visible with great detail. Based on the ex- 
perimental results the specific excitation of the two lobes 
could not be explained. However, the SNM amplitude 
maps as well as the spectroscopy of the duty cycle^^ gave 
a strong hint for possible in-plane rotational motion by 
a small angle. That was the starting point for our DFT 
studies. 
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FIG. 1: (Color online). Minimum-energy conformation 
of a copper phthalocyanine (CuPc) molecule adsorbed on 
Cu(lll). The molecular Cu atom is located on a bridge site 
and two of the molecular lobes lie along the dense atomic row 
in [-110] or equivalent direction (dashed white line) on the 
surface in agreement with the experimental findings^^. 



III. DENSITY FUNCTIONAL THEORY 
DESCRIPTION 

In order to understand the induced rotation of CuPc 
on Cu (111), DFT studies were performed to obtain a 
maximum of information about the electronic and geo- 
metric properties of the molecule on the surface. The 
calculations were done using the vasp code'^'^. 

The calculations were performed for a 9x10 Cu-atom 
unit cell with 4 layers. This large unit cell is required in 
order to reproduce dilute molecular densities. The large 
cell allows to use a single k-point. The PAW scheme 
for the atomic potentials^ was used, and the planewave 
basis set was expanded up to a cutoff energy of 300 eV. 

We evaluated the total energy of different high- 
symmetry conformations of the adsorbed molecule. 
Among all of them the minimum energy corresponds to 
the molecular Cu atom sitting on a surface bridge site 
and one of the molecular symmetry axis aligned along the 
[-110] direction (dense atomic row on the surface). Fig- 
ure[l]shows the minimum energy configuration, in p erfect 
agreement with the experimental observation&ISl. This 
adsorption conformation has also been found for CoPc 
on Cu(lll) in a recent joint experimental and theoreti- 
cal studySS as well as in other theoretical studie^^aiMI 



A. Adiabatic potential energy surface 

The PES as a function of the rotation angle, 0, between 
one of the molecular in-plane axes and the surface [-110] 
direction, was computed by fixing one of the pyrrole-N 
atoms and the molecular Cu atom at their relaxed po- 
sitions and rotating the corresponding interatomic axis 
to the desired angle. The structure was relaxed until all 
forces within the molecule and the first two substrate lay- 
ers fell below 0.02eV/A, keeping these two atoms fixed. 
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FIG. 2: (Color online). Adiabatic potential energy surface 
(PES) of a CuPc on Cu (111) as the molecule rotates around 
a surface normal defined by the position of its central Cu 
atom. Symbols correspond to computed values and lines are 
two different fits of the computed data. The two fits yield a 
frustrated rotation barrier between the equilibrium points at 
<j!> = 0° and <j) = 6.9° of 50 meV and 70 meV, respectively. 
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DFT+D2 
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2.574 


2.805 


E, 


hem (eV) 


4.21 


5.08 



TABLE I: Computed height distances z between the Cu 
central atom of the molecule and the surface plane and the 
molecular chemisorption energy Echem obtained with LDA 
and DFT-D2. 



phthalocyanines on noble metal surfaces^. In order to 
assess the accuracy of the LDA PES, we repeated the 
calculations just for the (p — 0°-case using DFT-D2 cal- 
culations including van der Waals interactiorPSl as imple- 
mented in the vasp code. Table ll] compares the results. 
The adsorption energy difference is 16%, and the differ- 
ence in adsorption distances is 8%. As is well-known, 
LDA compares favorably in these values with more real- 
istic methods. And, particularly in the present case, the 
LDA PES should be a good estimation of the adiabatic 
PES because the missing van der Waals force is a long- 
range interaction which is little affected by atomic details 
such as a small-angle in-plane rotation of the molecular 
axis. 



Afterwards, the constrain on the pyrrole-N was released 
and placed on an aza-N atom, repeating the ionic conver- 
gence to the same thresholds. In this way, the internal 
molecular and surface structure were relaxed for an ex- 
tensive set of ionic relaxation calculations, excluding the 
back-relaxation to the original equilibrium state. De- 
spite a careful convergence, we cannot rule out a small 
{A(j) « ±0.1°) uncertainty in the final positioning of the 
molecular axis. The PES shown in Fig. [2] reveals the ex- 
istence of a metastable adsorption position of the CuPc 
rotated by ±7° from the equilibrium position. This sec- 
ondary minimum is reminiscent of the geometry of full 
CuPc monolayers on Cu(lll) where the molecule appears 
rotated by ±7° from the substrate axis.'^^ESl 

The calculation of the rotational barrier height is of 
particular importance for the determination of the rota- 
tional rates as shown below. Figure [2] shows the adiabatic 
PES calculated with the above prescriptions. In order to 
smoothen the small calculation uncertainties a continu- 
ous function is fitted to the PES. Despite our efforts, we 
could fit different types of curves with different barriers. 
In the present calculations we have used two different 
curves corresponding to rotational barriers of 50 meV 
and 70 meV in order to study how critical the chosen 
PES is in the switching rate determination. 

All of these calculations have been performed within 
the local density approximation (LDA) . The rationale be- 
hind this is to use the overbinding error of LDA to ob- 
tain a physical molecule-surface distance compensating 
to some extend the neglected van der Waals interaction 
in LDA. It is well known, that van der Waals is a large 
component in the binding of large organic molecules like 



B. Electronic structure 

Transition metal phthalocyanines capture charge from 
noble metal surfaces^^. Similarly, in the case of CuPc, 
a full electron is captured into the first empty orbital of 
TT character (e2g in D4h notation). This has clearly been 
seen for CuPc on Ag(lOO)^. Copper surfaces are more 
reactive than silver surfaces, hence more charge is cap- 
tured. Our DFT results indicate a shift of the two e2g 
orbitals below the Fermi energy, and the overall trans- 
ferred charge approaches two electrons. Although the 
DFT Kohn-Sham orbitals strictly do not possess phys- 
ical meaning, they are commonly used to evaluate the 
molecule's charge population. The results can be under- 
stood as a qualitative prediction that can be easily put 
to test e.g. when computed constant current STM im- 
ages are compared to the experiment. Indeed, as shown 
in Ref. [15] there is qualitative agreement between the 
Tersoff-Hamann simulated image—' and the experimen- 
tal data. Here, we present calculations done within the 
Tersoff-Hamann theory but integrating over a given en- 
ergy window: 



/(x5]|V.(rO|Mi?-e.), 



(1) 



where E is the energy window from the Fermi energy, 
"01/ (r) is an eigenstate of the full system's Hamiltonian 
with its eigenvalue e,^. The window function F(E — e^) 
reads 



F{E-e,) = 



otherwise. 
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Figure [S] shows the evaluated constant current image 
using Eq. M for various values of E. For energy windows 
(equivalent to the applied STM bias) between -800 meV 
and 800 nieV,the resulting images slightly vary. The im- 
ages are dominated by the contribution of the extended 
TT-like e2g orbitals. Indeed, the density of states pro- 
jected on these orbitals-- show broad features spanning 
the above energy range. Due to a lack of wave function 
amplitude on-top of the molecular Cu atom, the STM 
image shows a depression at the molecular center. The 
characteristic four lobes resembling the two e2g orbitals 
appear as protrusions. 

Figure [3] also showcases a comparison between the cal- 
culated images for the molecule aligned with the [-110] 
axis (the ground state or minimum energy configuration) 
and the molecule rotated by ^ 7° (metastable state). 
Only small variations can be seen. There are some slight 
asymmetries due to inaccuracies during the coordinate 
relaxation, but more significantly the rotated molecules 
show sharper images (smaller sized spots in the lobe re- 
gion), indicative of a lower hybridization with the sub- 
strate. The simulated images also show a difference be- 
tween the lobes aligned with the [-110] axis and the per- 
pendicular ones. The latter show a larger local density 
of states (see for example the ground state image at 200 
meV, and also at -200 meV). This is due to the stronger 
interaction of the molecule with the surface along the 
[-110] axis, and has connections with the difference be- 
tween excitation probabilities as revealed in our excita- 
tion rate study below. Indeed, this finding relates to 
the experimentally observed factor of nearly two, sepa- 
rating the high-current and the low-current level of the 
RTN tunneling current signal. The experimental factor 
of about two was reproduced by our DFT study when 
comparing cross-sections within the calculated constant- 
height images. 



IV. TUNNELING-ELECTRON INDUCED 
ROTATION 

Energy transfer from a tunneling electron to an a dsor- 
bate has been studied in detail for various systems.'^SESl 
In the case of vibrational excitation by electron colli- 
sion, the resonant process associated with the transient 
capture of the electron by the target molecule has early 
been recognized to be very efficient for transferring en- 
ergy from an electron to a molecule. Resonant processes 
are very active in electron collisions on free molecules (see 
e.g.Hil) or adsorbed molecules^^, as well as in tunneling 
conditiona^^EH. For vibrational excitation the large mass 
ratio between electron and atoms makes a recoil mecha- 
nism only weakly efficient and any trapping phenomenon 
(a resonance) which significantly increases the electron- 
molecule interaction time can boost the vibrational exci- 
tation efhciency. The situation is completely different in 
the case of rotational excitation. The rotational motion is 
slow, so that an electron-molecule interaction can be seen 



as a sudden process with the fast electron colliding with 
a fixed molecule. The excitation process is then related 
to a recoil phenomenon in which the scattered electron 
transfers recoil angular momentum to the target. The 
specificity of the rotational excitation compared to the 
vibrational excitation comes from the change in the rel- 
ative orders of magnitude of electronic and nuclear mo- 
menta when going from linear to angular momenta. In 
the case of a free molecule, quantization of the angular 
momentum renders the angular momentum of a scattered 
electron at the same order of magnitude as that of the 
low rotational levels of the molecule, so that exchange 
can be very efficient and lead to significant rotational ex- 
citatiorpSl. In the present case, the non-spherical symme- 
try of the tunneling electrons and the frustrated rotation 
implies that neither the electron nor the molecule has a 
well-defined angular momentum. However, distributions 
of angular momentum can be associated with such sys- 
tems, and the efficiencyof the recoil in rotational excita- 
tion remains very higlPS. One can stress that magnetic 
excitation by tunneling electrons also involves angular 
momentum transfer. Essentially, spin excitation and ro- 
tational excitation can be described along the same lines 
and both can be very efficient^^ ^^. 



A. Scattering theory 

Let us consider the transition of an electron from an 
initial state \tpi) of energy e^ in the STM tip into a final 
state in the sample {ipf) of energy ef. The process may 
be inelastic such that e^ ^ ef. The excess energy may 
then be transferred to a frustrated rotational state of 
the molecule, defined as an eigenstate of the potential 
energy curve in Figure [2J The molecule is excited from 
its frustrated rotation ground state \Ro) to an excited 
state |i?n)j gaining an energy E — _E„ — Eq from the 
electron (£^0 and En are the energies of the initial and 
final rotational states). 

We can then say that the global system is initially in 
the state \i) = |i?o) ^ IV'i) at energy Ei = Eq + e^ and 
the final state and energy are [/) = |i?„) (E) \ipf) and 
Ef — En + ef , respectively. 

The transition rate per unit time for an inelastic elec- 
tron transfer from tip to the substrate is given by the 
T-matrix: 



1 



27r Y^ 



mj\^S{E,-Ef 



(2) 



In general, it is difficult to determine the connecting po- 
tential, and hence the T-matrix, in particular its inelastic 
part. However, in the present case, the molecular rota- 
tion is a slow motion compared to the electron trans- 
mission from the tip into the substrate. Therefore, one 
can use a sudden approximation in which the T-matrix is 
evaluated for fixed positions of the molecule with respect 
to the substrate (fixed 0-angle) and then used to com- 
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FIG. 3: (Color online). Constant current simulated images for different energy windows, Eq. ([T|. Both, the molecule 
aligned with the [-110] surface direction and the molecule rotated by ~ 7° are shown. An image variation due to the different 
hybridization of the molecule with the substrate is apparent between the two types of molecules. 



pute the T-matrix element between initial and final ro- 
tational states. This approach has been introduced very 
early in the treatment of rotational excitation of electron- 
free molecule collisionP^ and later also used for adsorbed 
molecules^S. Expression ^ can be further simplified in 
the TersofF-Hamann approximatiorpSl the T-matrix for a 
fixed position of the molecule is taken proportional to the 
electronic wavefunction evaluated at the tip apex (ro): 



T,j = (V',|f|^/)«CV'/(</',ro), 



(3) 



where i/'/((?!), tq) is the molecule-|-substrate wavefunction. 
It is computed for an STM tip located at fo, a fixed 
position with respect to the substrate and for an angle 
of the molecule rotated from its equilibrium position. C 
is the proportionality constant. We explicitly write the 
dependence of the wavefunction on the molecular angle (f> 
(relative position of the molecule and substrate) , for its 
later use in the evaluation of the rotational excitation. 

We can now replace each quantity in expression (pi) by 



its value 
1 



" Y^ ^ /t»p(e»)[l-/s«b(e/)] 

ij n>NR 
X \{Ro\{^^\f\^f)\Rn)\^S{Eo 



e,-E^- e/).(4) 



In the summation we explicitly defined electron flow 
from the tip to the substrate by using the Fermi factors of 
the tip flip and the substrate fsub- The index Nfi is that 
of the i?„ level above which the molecule actually rotates. 
Indeed, excitation from the ground state Rq to low lying 
rotational excited states i?„ is not sufficient to lift the 
molecule out of its equilibrium potential well (Fig. I2|. 
The rotational excitation must be strong enough for the 
molecule to overcome the potential barrier separating the 
equilibrium well from the two metastable ±7° side wells 
(see Fig. [2|. This limit corresponds to the level N^ in 
the formulas. For the PES with a rotational barrier of 70 
meV, Nfi^ is equal to 25 . This emphasizes, that jumps 
from one well to the other have to involve highly excited 
states of the frustrated rotation and thus, a very efficient 
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FIG. 4: Probability of de-excitation into the right weU from 
the level n of the rotational PES, Eq. ([8|. For levels below 
the 70 meV barrier there is an alternation of levels among left 
and right and central wells. Above the barrier, the rotational 
states become increasingly delocalized about all wells and the 
probability to find the molecule in the right well approaches 
1/3. 



rotational excitation process is needed. 

Making use of i) the Tcrsoff-Hamann approximation, 
ii) Eq. ([3| and Hi) the zero Kelvin expression for the 
Fermi distribution functions Q{ep — e), we can further 
simplify the expression by assuming that the tip is made 
of a material where the details of the band structure 
can be summarized in a density of states (DOS) func- 
tion Dtip{e). Finally, the rate for the excitation over the 
rotational barrier 1/Tj„e, is 
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(5) 



The total electron current is equal to the electron tran- 
sition rate from tip to substrate times the electron charge: 
I = e X -^—, where —^ is similar to Eq. dsl) except that 

TTot I Trot . . , O 

the summation over the rotation index n runs over the 



entire rotational spectrum. In the case of the total cur- 
rent, there is no restriction in the sum over n, while in the 
inelastic rate, Eq. (Is]), the summation is restricted to the 
over-barrier levels n > N^. We can define the fraction of 
the current that induces rotation from one potential well 
to another, also known as the inelastic electron fraction, 
by the ratio of the inelastic current to the total current: 



Tine TTot 



(6) 



The benefit from the implementation of the inelastic frac- 
tion is that the unknown constant C in Eq. ([3| and the 
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FIG. 5: General excitation and de-excitation scheme in be- 
tween rotational levels within the potential energy surface 
leading to the molecular frustrated rotation. A direct de- 
excitation from the highly excited level (upper red solid line) 
to the left side well is indicated by an arrow. The de-excitation 
can also take place via many intermediate states depending 
on the degree of excitation. Here we choose one intermedi- 
ate state (dotted line) to exemplify the mediated two-step 
de-excitation. 



tip's DOS factor cancel out under the assumption of a 
roughly constant DOS of the tip over the energy range 
of interest. 

Using the above expressions, the inelastic electron frac- 
tion becomes: 
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(7) 



The inelastic electron fraction ?/ does not evaluate the 
number of molecules that actually do rotate, but rather 
the fraction of electrons that excite the molecule over a 



certain rotational level Npi. Indeed, on the time scale of 
the experiment, a molecule excited by a tunneling elec- 
tron from the rotational ground state Rq to a rotational 



level Rn will quickly relax to a lower state, i.e. it will 
end up localized in one of the three potential wells and 
the experimental transition rate concerns transitions be- 
tween these potential wells. Relaxation from level i?„ can 
involve collisions with substrate electrons (electron-hole 
pair creation in a process very similar to the one dis- 
cussed here) or the transfer of the rotational excitation 
(coordinate cj)) to other degrees of freedom of the heavy 
particles in the molecule-substrate system (intramolecu- 
lar relaxation or transfer to phonons) . Figure p^ shows a 
simplified scheme of the excitation/de-excitation process 
leading to the transfer from the equilibrium potential well 
to another. For simplicity only one intermediate level in 
the relaxation is displayed in the figure. 

We did not try to evaluate in detail this very com- 
plex relaxation process but resort to a geometrical sta- 
tistical approximation. We evaluate the branching ratio 
for the de-excitation from the level Rn towards a given 
metastable well in terms of the weight of the rotational 
wavefunction i?„ in the 0-range of the metastable well, 
approximated as the 0-region beyond a critical angle </>£. 
By inspection of the PES, we choose 0^ = 5°. The frac- 
tion of molecules excited to a given level n that are even- 
tually trapped in the right well is: 



PR{n) 



15° 



\Rnm^dcf>/ 



15° 



and we can assume that the same probability rules the 
trapping in the left well. The ±15°-limits reflect that the 
molecule is not freely rotating above the surface, but it is 
confined to a 30°-sector due to the molecule-|-substrate 
joint symmetry and to the high barrier separating the 
equivalent regions (see Fig. [2]). Figure [4] shows the 
probability of de-excitation into the right well from the 
level n of the rotational PES, pii{n). For low n, the ex- 
cited states are localized in the equilibrium well so that 
the branching ratio is equal to zero. Above the energy of 
the bottom of the metastable wells, some of the states are 
localized in the side wells. Owing to our approach which 
considers even and odd functions, the branching ratio is 
equal to 0.5 for these states. Above the frustrated ro- 
tation barrier (70 meV for the shown case), i.e. for the 
states included in the summation in Eq. ([7]), the proba- 
bility steadily increases to 1/3, consistent with a roughly 
equal partitioning among the wells. 



|i?„(0)pd0, (8) 



The inelastic fraction of electrons that make the 
molecule jump from the equilibrium well into the right 
well, ?7fl, is then obtained by including the branching ra- 
tio PR{n) into the expression ([7]): 
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(9) 



At this point, one can stress that the quantity rjji can 
be directly compared with the equivalent quantity in the 
experiment, Ref. 1151 Furthermore, the excitation proba- 
bility r]ji{fQ, V) depends on where the tip apex is located, 
To, i-e. expression (|9| provides a geometrical map of the 
excitation process. 

As was experimentally showrPSl^ the derivative of the 
efficiency with respect to the bias gives density of states 
information. This can easily be comprehended, as the 
derivative changes one Q function in expression ([9]) into 
a delta function. The obtained quantity strongly resem- 
bles the local density of states with information on the 
involved rotational states. 

The above treatment parallels that used to compute 
the changes in conductance due to magnetic excitations 
(see a review in [27| . Indeed, both cases deal with an 
angular momentum transfer, in the present case, the ro- 
tational angular momentum of an adsorbed molecule, in 
the magnetic case the spin of a magnetic molecule. The 
theoretical treatments just evaluate the sharing of the 
electron flux among the different rotational or magnetic 
channels, respectively. 



B. Results 



Figure [6] shows the spatial distribution of the inelastic 
fraction of electrons producing the molecular rotation. 
The maximum corresponds to ~ 2.5 x 10~^ for a bias of 
0.5 V. The equivalent experimental numbepISl is rather 
~ 1.2 X 10~^ at 0.6 V. The calculations have been per- 
formed for a rotational barrier of 70 meV. As we show 
in Fig. [2j there is some uncertainty in the determination 
of the rotational barrier and 50 meV is also consistent 
with our computed PES. As expected, a change in the 
barrier height strongly influences the rate value. Indeed, 
the rotational rate, for a tip located at the maximum 
efficiency spot, changes by a factor of ^ 3 between the 
70 and the 50-meV barriers. This is probably one of 
the largest sources of uncertainty in the present calcu- 
lation. Another uncertainty is hidden inside the PES 
calculations of section IIA. It corresponds to the PES as 
a function of the angle </> when all the other coordinates 
have been relaxed, so that it includes an implicit adia- 
batic assumption for these coordinates, justified by the 
slowness of the rotational motion. In addition, in the 
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FIG. 6: Inelastic fraction of electrons inducing the molecular 
rotation. The probabilities of all states above the rotational 
barrier are added and a factor 1/3 is used for the branching 
ratio into the right well. 



presence of distortions of the molecule geometry as the 
angle (f) is varied, a more complex Schroedinger equation 
should be used than the pure rotation one used above to 
define the R„ wave-functions. In the present case, the 
distortion of the molecular frame are below 1 pm during 
the rotation, hence justifying the use of a single variable 
in the Schroedinger equation. 

The calculation involves the integration over angle, (p, 
the PES rotational levels and the electronic wavefunc- 
tion, ipf{(/),rQ). Since each angular value of ij}f{4>,ri^) in- 
volves a full self-consistent calculation of the molecule on 
the surface for that angle, it is a very costly calculation 
to perform. However, high i?„ rotational states exhibit a 
large number of nodes and thus require a large number of 
0- values for the integrals in Eq. ([9|. We balanced accu- 
racy and computation costs by discretizing the integra- 
tion with nine </)-values at which full DFT -calculations 
were performed. We then used a linear interpolation of 
ipf between the ah initio points. In this way we can have 
a double scale of integration points for the V'/ and i?„ 
parts of the integral. 

The calculation shows that the two molecular lobes 
aligned perpendicular to the dense atomic row direction 
[-110] yield the largest rotational rate, in excellent agree- 
ment with the experimenlflSl. The ratio of rates between 
the lobes perpendicular to [-110] and the parallel ones 
is roughly 5, Fig. [6] in very good agreement with the 
experimental one. By studying the electronic states of 
molecule and surface and their relaxed geometry, we con- 
clude that it is the local binding of the molecule to the 
surface that determines the efficiency of the angular mo- 
mentum transfer, i.e. the excitation mechanism is a local 
electronic effect. 



The DFT study on CuPc on Cu(lll) yields a molecu- 
lar adsorption conformation in good agreement with the 
experimental datcP^. The molecule is adsorbed flat on 
the surface with the central Cu atom sitting on a highly 
symmetric bridge site with two out of the four molecu- 
lar lobes aligned along the densely packed [-110], [0-11] 
or [10-1] directions. No buckled or bent conformation 
is found. In addition, we calculated the evolution of 
the total energy of the system when the molecule is ro- 
tated around a surface normal, going through the central 
atom of the molecule. The resulting PES reveals the ex- 
istence of two local minima for symmetrical rotations by 
±7° clockwise and anti-clockwise from the equilibrium 
position. The molecule can undergo a frustrated rota- 
tional motion around its equilibrium position as soon as 
it gains an energy higher than the rotational barrier (70 
nieV) by flipping from the equilibrium well into one of 
the side wells of the PES as depicted in Fig. [2] These 
computational results confirm the experimental interpre- 
tation that the blurred STM images of one opposite pair 
of benzopyrrole rings of CuPc are caused by jumps of 
the molecule between three adsorption wells induced by 
tunneling electrons^. 

The simulated STM images show that the four CuPc 
lobes are not equivalent but split into two groups (lobes 
along the dense atomic rows and perpendicular to it). 
Their position relative to the underlying Cu lattice are 
different and consequently they hybridize differently with 
the substrate. The strong hybridization of the 7r-like or- 
bitals located on the molecular lobes leads to a sizable 
charge transfer of 2 electrons from the surface to the 
molecule. As a result, the simulated STM image is dom- 
inated by the contribution of the e2g (7r-like) orbitals of 
the molecule. This contribution is slightly different on 
the two types of lobes in good agreement with the sym- 
metry reduction observed in the experiment. 

By analyzing the telegraph noise in the tunneling cur- 
rent, which causes a fuzzy appearance of CuPc in STM 
images, Ref. [TS] showed that the molecular rotation rate 
is maximal on the two lobes perpendicular to the dense 
atomic row [-110] direction, whereas rotation cannot ef- 
fectively be induced on the remaining lobes. We have 
performed calculations that revealed the excitation mech- 
anism of tunneling electrons inducing the frustrated rota- 
tional dynamics and leading to fluctuations in the STM 
tunneling current signals. Our approach is based on the 
sudden approximation for the rotation and the Tersoff- 
Hamann approximation for the T-matrix in the tunneling 
regime. The main result of the theory is that the tunnel- 
ing electron has a finite probability of exciting the molec- 
ular rotation: the electron is briefly in contact with the 
molecule and departs leaving it in one of the rotational 
levels; when the excited rotational level is above the po- 
tential energy barrier separating the three frustrated ro- 
tation wells, then the jump of the molecule from one well 
to another becomes possible. This process is analogous 



to the one found in magnetic excitation by tunneling elec- 
trons, where the adsorbate spin changes depending on the 
we ight o f each magnetic state in the tunneling symme- 
tryESET] Q^j. calculations show a clear difference in the 
excitation rates over the two types of molecular lobes, 
identical to the selectivity observed experimentally. We 
connect this selectivity with the slight differences in the 
STM images of the molecular lobes, attributed to the 
different hybridization between the e2g orbitals (mainly 
located on the molecular lobes) and the substrate. These 
different hybridizations lead to different rotational rates. 
The lobe selectivity observed in the telegraph noise in- 
duced by tunneling electrons thus appears to be of elec- 



tronic origin. 
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